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COMPUTATIONAL FLUID DYNAMICS AND AEROTHERMODYNAMI CS 
(Development o-f Nonequilibrium Models Applicable to AOTV/AFE and 

Super-orbital Flight Regimes) 

I . Introduct ion 

This report covers the period -from January 7, 198? thru June 30, 
1987. It will brie-fly decribe and summarize the current status o-f the 
project, the research e-f-forts during the last six months, and projected 

-future e-f-forts. In addition, copies o-f viewgraphs recently used during 

I 

a brie-fing to Johnson Space Center (JSC) personnel are included as an 

l 

appendix. 

1 1 . Status and Personnel 

Officially, the present project was awarded on January 7, 1987. 
However, final paperwork was not completed until the middle of February; 
and by this time personnel assignments for spring semester had already 
been decided. Consequently, no graduate students were assigned to this 
project until April 1, 1987. Due to this delayed start on the project, 
it may be necessary to request a no-cost extension to the project 
sometime in the next few months. 



At the present time, the staff associated with the project are: 
Leland A. Carlson — Principal Investigator 

Approximately 1/4 time 

Glenn Bobsk i 1 1 — Graduate Research Assistant 

Apr i 1 - May, Hal f t ime 
June to present, Full time 
Thomas Gaily — Graduate Research Assistant 
June to present, Half time 
Robert Greendyke — Graduate Research Assistant 

June to present, Full time 

The research work associated with this project will -form the basis 
■for the Masters thesis work o-f Mr. Bobsk i 1 1 and Mr. Greendyke. Mr. 
Gaily received his Masters degree in May 1987 and is currently working 
on his doctoral degree. His current research efforts are split between 
this project and one involving transonic wing design and analysis. 

111. Research Prooress 

At the present time the e-f-forts on this project are organized into 
two areas. The first involves Tom Gaily and Dr. Carlson and is directed 
towards the development of an approximate stagnation point solution. 
The second is directed towards approximate flowfield studies which can 
be used to develop and investigate shock jump, electron temperture, 
radiation, vibration-dissociation coupling, and chemistry models. The 
individuals associated with this second portion of the project are Dr. 
Carlson, Bob Greendyke, and Glenn Bobsk i 1 1 . 
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Approximate Stagnation Point Solution Studies 

The objective of this effort is to develp a rapid and "accurate" 
approximate stagnation point solution which can be used for parameteric 
and/or design studies and which includes approximately the effects of 
chemical and thermal nonequ i 1 i bi rum, thermal conduction, radiative 
gasdynamic coupling, and convective and radiative heat transfer. Since 
this model is currently under development and undergoing daily revision, 
pressentation of details will be delayed until later reports. However, 
a general outline of the approach is presented in the appendix. At the 
present time, computer solutions complete with graphics are being 
obtained in an interactive fashion; and the entire approach appears 
promising. 

An example of the type of results which can be obtained with the 
present approximate solution is shown on Figure 1. This case is for a 
freestream velocity of 1<4 km/sec, an altitude of approximately 63.5 km., 
and a shock detachment distance of 5.88 cm. As can be seen, the flow 
rapidly equilibrates behind the shock wave and is in chemical 
equilibrium over most of the shock layer. Also, the flow exhibits some 
radiative cooling and coupling as evidenced by the steady decrease in 
temperature and degree of ionization in the central equilibirum part of 
the shock layer. 

Now, normally one would expect that radiative heating from the 

* 

transparent portion o-f the spectrum would correlate with the -freestream 
parameter 0 .5*<Freestream Densi ty)*(Freestream Uelocity Cubed), which is 
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a measure of the kinetic energy -flux or power per unit area in the 
freestrearo, and scale with the shock detachment distance. 
Interestingly, this case, under that logic, corresponds closely to the 
condition o-f peak heating observed in the FIRE II -flight experiment even 
though that case was at about 10.48 km/sec and 53.7? km. The present 
solution does contain an approximation -for the heating -from the visible 
portion o-f the specturm, which is essentially transparent; and when this 
value is adjusted to account for differences in shock stand-off 
distance, the result is about 350 watts/sq.cm. While the project has 
not yet obtained the actual FIRE numerical data, graphical plots of the 
measurements indicate that stagnation radiative heating values of 350 to 
400 watts/sq.cm, were measured from the visible wavelength regions at 
the time of peak heating. 

A result obtained from the current approximate method which 
corresponds more closely to a FIRE II trajectory point is shown on 
Figure 2. This result is in reasonable agreement with detailed 
computations recently presented by Bird (Reference 1) and obtained using 
the direct simulation Monte Carlo method. Again, the flow behind the 
shock front equilibrates relatively quickly and most of the flowfield is 
in chemical equilibrium. However, due to the lower freestream velocity, 
the degree of ionization is relatively low and there is very little 
radiative gasdynamic coupling. 

As mentioned previously, the present method contains *a highly 
approximate method of estimating the stagnation point heating due to 
radiative transfer from the visible region. J n addition, a version of 
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the code is being developed which uses the resultant -flow-field solution 
to compute a more accurate estimate o-f radiative heating. The latter 
approach uses a -five step radiation model which includes emission and 
absorption -from the vacuum ultraviolet <VUV) continuum and lines, 
visible continuum, and in-frared <IR) lines. Figure 3 compares the 
present predictions with data -from the FIRE II experiment. The 
measurements -from FIRE II were only in the visible plus IR spectrum and 
only those results are shown. The present results are plotted over a 
horizontal band since the present trajectory point did not match exactly 
a FIRE II point. Thus, one end o-f the spread corresponds to a match in 
freestream velocity while the other corresponds to a match o-f altitude. 
Interestingly, the present results "surround" the experimental data. 

However, in this investigators opinion, this good agreement should 
not be viewed as verification of the method but rather as an indication 
that it is relatively easy to match the FIRE data in the plotted 
wavelength region. Phrased differently, the radiative heat flux in the 
visible plus IR is not a sensitive indicator of the accuracy of a 
flowfield calculation or model. 

Nevertheless, it is believed that the present five-step radiative 
heat transfer results are int resting. These results, which should be 
viewed as very preliminary, are as follows: 
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Wav el ength 

Type 

Heat Transfer 

(Angstroms) 


(Watts/sq .cm.) 

620-1100 

VUV Continuum 

0.035 

1100-1300 

VUV Continuum 

51 .89 

1300-1570 

VUV Lines 

1219. 

1570-7870 

Visible Continuum 

18.08 

7870-9552 

IR Li nes 

108.5 


First, it should be noticed that in the visible plus IR region, 
most of the heating is from IR lines. This trend is in accord with 
shock tube data and calculations o-f other investigators. However, the 
most surprising feature, is the strong contribution from the VUV and in 
particular from VUV lines. While, the contribution from this region 
will depend strongly upon shock detachment distance <i.e. optical path 
length) and the degree of nonequilibrium, these results indicate that 
VUV radiation and not the visible may dominate. While this possibility 
has been ignored by many investigators (assuming that the VUV will all 
be absorbed and not reach the body), it is in agreement with the recent 
calculations of Bird and is implied in the comparison between the 
visible plus IR and total heating measured on FIRE II. 

Considering the uncertainty of these results and the lack of shock 
tube data in this wavelength regime, consideration probably should be 
given to instrumenting the AFE to either measure heating from the VUV or 
to provide measurements from which it can be deduced. 


During the next portion of the project, it is planned to continue 
development o-f the approximate stagnation point model. Tentatively, the 
■following tasks will be conducted: 

<1) Incorporate methods -for accurately computing radiative heat 
transfer using a -five step absorption coefficient model. These 
techniques will include the effects due to chemical and thermal 

nonequ i 1 i br i urn . 

(2) Develop and include in the method appropriate techniques for 
computing the covective heat transfer to the stagnation point. 

(3) Improve the chemistry correlations used in the method. 

<4) Automate the method so that it is easy to use and suitable for 
parametric type studies. 

(5) Consider and if possible include slip boundary conditions so 
as to extend the region of validity of the model. 

( 6 ) Use the model to perform parametric type studies and compare 
results to experimental and flight data. 

Approximate Flowfield Solution Studies 

The primary objective of this portion of the program is to develop 
and study a variety of of models and approximate models for various 
physical phenomena associated with the flowfield. In order to develop 
such models, delineate their applicability, and determine the 

consequences of their usage, it is essential to have a rapid flowfield 

* 

solver. After considerable investigation, it was decided to use for the 
flow solver a chemical and thermal nonequilibrium Maslen type of inverse 
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solver due to Grose (Reference 2). In this method, a shock shape is 
assumed and the -flow-field and body shape are determined as part o-f the 
sol ut i on . 


While a Maslen model utilizes many assumptions, experience by this 
and other investigators has shown that such models yield very good 
engineering type o-f results. In addition, this model, which is based on 
the Euler equations, has several distinct advantages to the present 
project. First, it is computationally very -fast. Second, it can be 
utilized to compute the solutions along only a -few selected streamlines 
•for those situations in which an entire -flow-field solutions is not 
needed. This aspect is very important currently since many studies will 
be needed to determine the e-f-fects etc. o-f using various approximate 
models. Finally, the selected Maslen approach should be able to model 
AOTV/AFE type o-f vehicles reasonabl y wel 1 since they are basically blunt 
sphere-cones. 

Currently, the model has been coded and test runs have been 
conducted. These test cases have included chemical nonequilibrium 
< seven species and six reactions) as well as vibrational nonequilibrium, 
and some typical results are included in the Appendix. It is planned to 
use this -flow solver to study shock jump, electron temperature, 
radiation, vibration-dissociation coupling, and reaction chemistry 
mode 1 s . 

* 

One o-f the primary problems in solving chemical nonequ i 1 br i um blunt 
body -flows with the Euler and/or Navier Stokes equations is the 
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disparite time scales between various phenomena. Typically, some o-f the 
chemical reactions have characteristic times much shorter than others 
and much shorter than the -flow characteristic time. These differences 
make the resultant equations stiff, which often leads to stability 
problems and long computing times. However, from a physical standpoint, 
chemical reactions in air often proceed in steps. For example, oxygen 
dissociation occurs much more rapidly than nitrogen dissociation, and 
both occur faster than ionization. Thus, in various flight regimes it 
might be possible to essentially uncouple the processes by assuming that 
one or more of them is complete in the "shock front" and incorporting 
the effect into the shock jump conditions. 

Consequently, work is currently in progress under this portion of 

the project to develop several shock jump approximations, incorporate 

them into the flowfield solver, and study their applicability, 

subsequent effects on the flowfield, and the computational effort <or 

savings) involved. Currently three types of shock jump conditions, in 

addition to the default condition of frozen mass fractions, are being 

developed. The first, under the assumption that conditions are such 

that oxygen dissociation occurs essentially immediately behind the 

front, will determine post shock conditions assuming that 02 and 0 are 

in equilibrium with each other and that N2 and all other species 

concentrations are frozen. It is believed that this approximation will 

be valid over the blunt body portion of much of the AQTt) flight regime. 

* 

The second one is applicable only to pure N2 and assumes that 
nitrogen dissociation occurs rapidly and that N2 and N are in 
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equilibrium immediately behind the front. While not directly applicable 
to ATOU atmospheric flight, the use and study of this approximation will 
enable some rapid studies of superorbital velocity flight and should 
give insight into a variety of phenomena which are dominated or strongly 
influenced by ionization and nitrogen atoms. At speeds above 10 km/sec 
nitrogen atoms and ions are the primary species affecting radiative 
phenomena . 

Finally, a shock jump approximation is being developed which 
assumes that both nitrogen and oxygen dissociation are essentially 
complete imemdiately behind the shock front and that N2 and 02 are in 
equilibrium with N and 0. Again the validity and applicability of this 
approximation will be investigated. 

Another area which must be modeled appropriately in order to 
accurately predict AOTWsuperor i btal flowfields is the electron 
temperature, which is the temperature which controls radiative emission 
and absorption. At velocities above 10 km/sec there are sufficient 
number of free electrons that theoretically the electron temperature 
should be determined by solving an electron energy equation. However, 
the inclusion of this highly coupled equation into the problem greatly 
increases the computational time required fc' solution. Consequently, 
an investigation and evaluation of various methods of approximating the 

electron temperature for this flight regime is being conducted.. 

* 

On the other hand, for velocities be lew approximately 10 km/sec, it 
is generally accepted that the electron temperature is closely coupled 



to the nitrogen vibrational temperature; and some researchers have 
suggested approximating the electron temperature by setting it equal to 
the nitrogen vibrational temperature. While this approach is quite good 
over the latter portion of the chemical relaxation zone behind the shock 
wave, previous experimental and theoretical studies by this author 
indicate that the peak electron temperature may be si gn i 4 i cantl y lower 
than the peak nitrogen vibrational temperature and that the time of the 
peak may be later. These differences could be important in computing 
the resultant radiative emission. Consequently, various electron 
temperature models are being studied and will be incorporated into the 
flowfield program and evaluated. 

In conjunction with the evaluation of electron temperature models, 
various radiation models are also being studied. Since this effort has 
just been initiated, details are not yet available. However, for 
temperatures above 8000 to 10000 K, it appears that step models based 
upon nitrogen atmos and ions will be adequate. However, such models 
might possibly be simplified, and this possibility will be investigated. 

At lower velocities and/or temperatures molecular band radiation 
becomes important and current project studies have not definitely 
determined whether or not step models can be used. The work in th s 
area will, of course, continue. 

By far one of the most important areas being studied currently 
under the project is vibration di sssoc i at i on coupling models. Vibration 
dissociation coupling strongly influences species concentration 
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variations, vibrational temperatures, and hence electron temperatures 
and radiative emission-absorption. Currently, the flow-field program has 
a coupled vibration dissociation (CVD) model; and this approximation is 
being modified to a coupled vibration dissociation vibration (CVDV) 
model. Eventually, it is planned to include in the program the option 
of using either of these two models along with the CVDV preferential 
model and the more recently proposed Park model (Reference 3). In 
addition, present project studies indicate that the choice of 
vibrational relaxation data and correlations can have a significant 
effect on results. Consequently, in the next period the similarities, 
differences, and consequences of using each of these models along with 
various relaxation data will be investigated for the AOTV regime. 

In conjunction with vibration dissociation coupling efforts work 
has also started on investigating the effects of various reaction 
chemistry models. It is planned not only to perform parametric studies 
of various rates but also to determine the effect on the flowfield and 
heat transfer of various proposed reaction rate temperature dependence 
models. While historically reaction rates have used heavy particle 
translational temperatures, several investigators have proposed for 
various reactions using electron temperatures, vibrational temperatures, 
or combinations of various temperatures. Current studies indicate that 
these various choices can affect rates and the extent of nonequilibrium 
by orders of magnitude. Thus, work is in progress to modify the flow 
code to allow various temperature options. After this work is 
completed, studies will be conducted to determine the effect of various 
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temperature dependence models. It should be noted that -for these 
studies the -fact that the inverse -flow code can be used to compute only 
a -few streamnlines will result in significant time savings. 


IV. Future Efforts 

As indicated above and based upon discussions during a recent 
present to NASA personnel, work during the next project period will 
continue the development of the various approximate codes. During the 
next reporting period these will be used as follows: 

<1) To perform prametric studies of shock jump, electron 
temperature, and radiation models. 

<2) To perform parametric studies of vibration dissociation 
chemistry coupling and reaction chemistry models. 

<3> For the evaluation of various approximate models. 

<4) For the evaluation of various approximations associated with 
stagnation point solutions and to perform parametric studies of the 
influence of various parameters on stagnation point heating. 

Finally, it should be noted that the present effort is directed not 
only towards the AFE/AQTV velocity regime but also towards velocities 
above 10 km/sec. There are essentially two reasons for considering the 
higher velocities. First, investigation of superorbital phenomena and 

the development of engineering models for that regime will be of value 

* 

for future lunar and Mars missions. Second, the entry velocity of 10 
km/sec for the AOTV should be viewed only as a "maximum cruise” point. 
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As -for any aerospace vehicle, the actual -flight envelope used to 
determine the vehicle design must, in order to account -for emergencies 
and various mission options, encompass a larger velocity and altitude 
range. Thus, those portions o-f the current studies which are directed 
at velocities above 10 km/sec are -from an engineering standpoint also 
applicable to an AFE/AOTV type o-f vehicle. 
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CONPUTATI ON A l_ FL_U I O 
DYNAM I 0 3 AND 
o r M E RM O DYNAM I O S5 

(Nonequ i 7 i br i um Mode 1 s- Ap p 1 i cabl e to AOTV/AFE 

F7 i ght Regimes) 


PRO JE C~T STATUS 
E-f-fective Start Dates 1 April 87 


TW O F>"T S 


PART I : 

Approximate Stagnation Point Solution 
<Gally and Carlson) 


PART I I s 

Approximate Flow-field Solution Studies 

<a> Shock, Te , and Radiation 
Appr ox i ma t i on s 
<Greend/ke and Carlson) 

U i b . - Temp. Coupl i ng Models, and 
Effect on Chemistry 
(Bobsk i 1 1 and Carlson) 


<b> 


i---* f 3, f> r: ox i i - irvt* i:-_ str* gnat i on 

R* O I N'F SOLUT I OFJ 


OBJECTIVE : 

Develop an "accurate" approximate 
solution suitable -for parametric studies. 

STARTING POINTS 

Superorbital velocities and continuum 
-flow. Later extend to lower velocities and 
transition region. 

FLOWFI ELD STRUCTURE: 


1 1 

Con due t i ng 

I . Conduct i ng 

l\ 

1 1 


1 

l\ 

1 1 

Chem. Noneq. 

1 Chem. Eq . 

l\ 

1 1 


1 

l\ 

t i 

Complete Dissoc. 

1 Complete Dissoc 

l\ 

I i 

in Shock Front 

1 

l\ 

1 1 


1 

l\ 

1 1 

Ambi polar Diffusion 

1 No Mol . Diff . 

l\ 

1 1 

.1 

1 

IN 

1 i 

u = 0 

1 u = zero 

IN 

1 1 

p = constant 

1 p - constant 

IN 

1 i 


1 

IN 

1 1 

Radiative Losses 

1 Radiative Losses from 

IN 

1 1 

Dominated by UUU 

1 Optically Thin Part ofl\ 

1 1 


1 Spectrum 

IN 

1 1 


1 

IN 

1 1 

TH £ Te 

1 TH = Te 

IN 

1 1 


1 

IN 

1 1 

Te Constant 

1 

IN 


Bow Shock End of Noneq. Zone Wall 


GOVERNING EQUATIONS 



21 




CHEMISTRY s 

"Uncouple" Chemistry by Assuming 

P&j j) ■£ £ — £ ^ h , iy 5 

-for n o'a) 

£ = £ ecj Cl- exp < - B t r >3 

whore B is expt./empi rical 
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THICK EQUATION: 


A 

h " - 

RePr 


A = K, L 

't 


D — 


< 1 


*r = 


B 


'Ey > h" = 
A 

-z.tr 

e + C 

< 1 - *r 

) c 

A -*& 

Die 

T 

A 


RePr 

2 tt s r 4 

o 

II 

V 

- tv* > G- €/s • B 

* $ 

f>< u c h * 




Re = 

■ P u ■ h 
r 4 sc 


1 

II 

|C 


SHOCK B.C. s 

Includes- thermal conduction 

Actual value coupled to the solution 


THIN ZONE 


ABSORPTION COEFFICIENT 


K i “ K <. 


< ? 'P. 


-t 


RADIATIVE TRANSFER: 


d q 


dx 


K 


t 


f 

A 




where n is typically 3 



CHEMI STRY : 

Assume £ = £ < h > and use 1 inear 
s e Qrne nts. 
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WALL B . 


C . : 

Currently assuming “hot" noncatalyt ic 

cm a 1 1 . 


THIN EQUATION: 



PHYSICAL COORDINATES: . 

Obta i n © d f r om solution and 
trans-format i ons. 
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F l_l~T URE EL R F-" O F iT S 

1. Accurately compute radiative heat transfer 
using a 5 — step model . 

2. Compute covective heat transfer. 

3. Incorporate nonequilibrium effects on 
absorption coef f i c i ents. 

4. Improve relaxation “distance" values. 

5. Improve chemistry cor r e 1 a t i on s . 

<5. Improve match points and thick downstream 
boundary condition. 

7. Automate 

8. Compare to data and perform parametric 
s t ud i e s . 

9. Extend -further into AOTV/AFE -flight 
r e g i me . 


29 



> 


s-=. RF>R! ox I MAI' E 
S O l— l_*T I ON 


I- I ow \~ I El — L> 

ST l_l O I E S 


FLOWFIELD MODEL 

Us i ng a nonequ i 1 i br i um Masl en type o-f mode 1 
due to Grose. 

Advan t ages s 

1 . Fast 

2. Can do individual streaml ines 

3. Should model AOTU/AFE reasonably 
we 1 1 

Usage : 

1. Study Shock, TE, and Radiation 
Mode 1 s 

2. Study Vi brat i onal -Chemi stry 
Temperature Coupling models 

3. Study reaction chemistry models 


SHOCK JUMP APPROXIMATIONS 

Currently developing -four di-f-ferent cases: 

< 1 > Fr oz e n 

<2> 02 & O in equ i 1 N2 -frozen 

< 3 > Pure N2 with N2 & N in equil. 

< 4 > N2 , 02 in equil. with N & O 

Plan to determine where approx, appl icable, 
subsequent effects on -flow-field, and 
computational effort involved (savings) . 
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TEMPERATURE (1000* K) 


ELECTRON TEMPERATURE MODELS 


At high temperatures, can TE be 
approximated or must a TE DEQ be solved? How 
can problem be simpl i -f i ed? 

When significant number of molecules 
present, TE is coupled to TL*N2. But how is 
it coupled? Is TE®=T*JN2 okay? Is a more 
detailed model needed? 


Previous results — 



l’ic. 3. Comparison of theoretical prediction and experimental 
data for electron temperature. 
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RADIATION MODELS 


At high temperatures., above 8000 to 10000 
K, -flow consists of atoms and ions. Step 
models <5> appear adequate, but possibly can 
be simplified. Radiation gasdynamic coupling 
probably important. 

When molecular radiation important, can 
step models be used? How many steps? 
Radiation coupl i ng probably not important. 


UI BRAT I ON-CHEMISTRY COUPLING MODELS 

Var i ous possi bilities exists 

1 . CL>DL> 

2. CVDV Preferential 

3. Park Model 

Influence of vibrational relaxation data 
and cor r e 1 at i ons? 

The similarities, differences, and 
consequences of using each of these models in 
the AOT*v> regime needs to be determined. 


REACTION CHEMISTRY MODELS 

Must be studied in conjunction with 
vib-chem. models. 

Study influence of number of species and 

number of reactions on f 1 owf i e 1 d structure 
and computational effort. 

Perform parametric studies on rates. 

Study effect of various proposed reaction 

rate temperature dependence models. 

Significant information probably could be 

obtained by studying only two streaml ines. 
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UT L_» F-i EH O I F: ET OT I ON S 


1 . Parametric studies of shock jump, electron 
temperature, and radiation models. 

2. Parametric studies of vibration chemistry 
coupl i ng and reaction chemistry' models. 

3. Evaluation of various models. 

4. Continue to develop approximate stagnation 
point sol u t i on . 

NOTE: Are these acceptable? We want to 

provide results that are of value to you! 
Guidance Please! Duration? Renewal? etc. 
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